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Abstract 

Adopting the phantom (ghost) scalar field description of dark energy, 
we construct a general class of exact interior solutions describing mixed 
relativistic stars containing both ordinary matter and dark energy in 
different proportions. The exterior solution that continuously matches 
the interior solutions is also found. Exact solutions describing extremal 
configurations with zero ordinary matter pressure are also constructed. 
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The astronomical observations of the present Universe provide evidence for the existence 
of a mysterious kind of matter called dark energy which governs the expansion of the 
Universe [1] , [2] . The dark energy exhibits some unusual properties such as negative pressure 
to density ratio w (in hydrodynamical language) and violation of the energy conditions. The 
ratio w may be less than —1/3 corresponding to violation of the strong energy condition, 
or even less than —1 which is violation of the weak or null energy condition. Current 
experimental data shows that w is in the range —1.38 < w < —0.82 . 

The fundamental role that the dark energy plays in cosmology naturally makes us search 
for local astrophysical manifestation of it. In the present work we consider models of rel- 
ativistic stars containing not only ordinary matter but also dark energy. The existence 
of dark energy makes us expect that the present-day existing stars are a mixture of both 
ordinary matter and dark energy in different proportions. The study of such mixed objects 
is a new interesting problem and some steps in this direction have already been made (see 
for example [5]- [13] and references therein). 

A possible theoretical description of the dark energy is provided by scalar fields with 
negative kinetic energy, the so-called phantom (ghost) scalars [3], [1]. The negative kinetic 
energy, however, leads to severe quantum instabilities and this is a formidable challenge 
to the theory. However, there are claims that these instabilities can be avoided [14J. In 
general, the problem could be avoided if we consider the phantom scalars as an effective 
field theory resulting from some kind of fundamental theory with a positive energy [15], [IB] . 
In this case the phantom scalar description of the dark energy is physically acceptable. In 
this context it is worth noting that the phantom-type fields arise in string theories and 
supergravity [T9]- [22] . 

In this work we adopt a description of the dark energy by a phantom scalar. Then the 
Einstein equations in the presence of dark energy read 

Rnu = 87r(T^ - -Tg^) - 28^8^, (1) 
V M V>V = 4irp D . 

Here T^ v is the energy-momentum tensor of the ordinary matter in the perfect fluid descrip- 
tion with energy density p and pressure p: 

= (p + pju^Uu + pg^. (2) 

For the ordinary matter we impose the natural conditions p > and p > 0. The density 
of the dark energy sources is denoted by p D . The dark energy sources will be called dark 
charges. 

Since a very little is kown about the interaction of dark energy with the normal matter, in 
our model we consider only minimal interactions for the phantom field - we have not included 

*The perfect fluid description of the dark energy also suffers from instabilities due to the imaginary 
velocity of the sound. From a classical point of view the massless phantom field is even more stable than 
its usual counterpart [T7], [15] . 
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terms describing non-minimal interaction between the phantom field and the normal matter 
in the field equations dJJ). 

When we consider the local manifestation of dark energy on astrophysical scales, i.e. 
scales much smaller than the cosmological scales, the phantom potential U(<p) can be ne- 
glected and that is why we set U((p) = 0. This is not a principle restriction. Exact solutions 
with the desired properties as those presented below can be also found in the presence 
of phantom potential and non-minimal interaction between the ordinary matter and the 
phantom scalar. Since we do not seek mathematical generality but we are interested in the 
physics of the model we shall restrict ourselves to the case U((p) = and minimal interaction 
with the ordinary matter. 

In what follows we give exact solutions describing mixed relativistic dark energy stars. 
The dark star solutions are characterized by the mass M, the dark charge D and the 
(coordinate) radius R. 

We consider static, spherically symmetric and asymptotically flat space-time with a 
metric 

ds 2 = -e 2U dt 2 + e~ 2U+2X [e-^dr 2 + r 2 (d6 2 + sin 2 Odtf)} , (3) 

where the metric functions are functions of the radial coordinate only. The phantom field 
ip and the dark charge density p D are required to be static and to depend on the radial 
coordinate r only. For the normal matter fluid we impose the usual conditions for staticity 
and spherical symmetry p = p(r), p = p{r) and u^dx^ 1 = —e u dt. Applying the mathematical 
techniques developed in [23] we can generate exact interior solutions to the field equations 
(JT]) using any exact interior solution of the ordinary Einstein-perfect-fluid equations as a 
seed. More precisely, we have the following proposition: 

Proposition Let 

ds 2 E = -e 2X dt 2 + e~ 2x dr 2 + r 2 (d6 2 + sin 2 6d<J) 2 ), 

p E =PE(r), (4) 

Pe = Pe{t) 

be an interior solution to the ordinary Einstein-perfect-fluid equations, then the metric (dj] 
together with the functions 



e 2U = e 2Xcosh -P 

p = e 2U ~ 2X [p E cosh/3 + 3(cosh/3 - l)p E ] , 

P = e 2U ~ 2X PE, (5) 
Pd = e 2U - 2X {p E + 3p E ) sinh/3, 
ip = sinh j3\ 

form an interior solution to the field equations (QP where ft is an arbitrary real constant. 
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It is worth noting that the solution generating method described in the proposition can be 
applied for an arbitrary equation of state of the ordinary matter. 

If the seed solution has a well-defined boundary r = R, where by definition Pe{R) = 
0, the same is true for the solution (jSJ), namely p(R) = 0. Therefore R can be interpreted as 
the coordinate radius of the dark energy star. On the dark star surface r = R, our interior 
solution (jSJ) matches continuously the following exterior solution 



2m\ cosh/3 ,„ / 2m^ cosh/3 



t2 



dsi xt = - ( 1 ] df + ( 1 



,/r " +r 2 (^ 2 + sin 2 



2m 
r 



1 / 2m \ 

Vext = ~ sinn P m 1 • ( 6 ) 



2 V r J 

For /3 = this solution reduces to the vacuum Schwarzschild solution describing the exterior 
gravitational field of a static and spherically symmetric star with mass m. The exterior 
solution can be easily and elegantly generated from the vacuum Schwarzschild solution if 
the 50(1, 1) symmetry of the dimensionally reduced field equations (jTJ) with p = p = p D = 
is used. However, it should be mentioned that ([6]) was obtained for the first time in [23] 
(with another method) and can be considered as a phantom counterpart of Fisher's scalar 
vacuum solution for a normal scalar [25], [26J. 

The mass and the dark charge of the dark star are given by 

M = ~f R t tV ^d 3 x= [ ( p + 3p) v ^d 3 x= (7) 

47r J Star J Star 

cosh/3 / £ 2X (pe + 3pE)y/—g~E~d 3 x = cosh /3m, 

J Star 

D = [ p D V^gd 3 x = ^-l V^d^ = (8) 

J Star 4?r JSL 



sinh/3 / e 2X (p E + 3p E )y/—gEd 3 x = sinh/3m 

J Star 



The same expressions for the mass and the dark charge are also obtained from the asymptotic 
expansion of the exterior solution. The mass M and the dark charge D satisfy the inequality 
M > \D\ as one can see from the above expressions. The dark charge D is a measure of the 
content of dark energy in the star. When there is no dark energy in the star (i.e. D — 0) 
our solution reduces to the seed interior solution describing ordinary perfect fluid star in 
general relativity. The initial solution parameters m and /3 can be expressed in terms of M 
and D as follows 



m = VM 2 - D 2 , cosh = M (9) 



VM 2 - D 2 
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In the case of ordinary stars (D = 0), the radius of the star R and the gravitational 
Schwarzschild radius Rq = 2m satisfy the well-known Buchdahl inequality [27J 



2m 8 



R 



<9' 



(10) 



For the mixed dark stars we have 



2M 2m cosh 



R, 



•ph 



R 



2m \ 
~R) 



(cosh/3-l)/2 



(11) 



where R p h — R (l — cosh/3 - )//2 is the physical radius of the dark star. By examining the 

right hand side of eq. (fTT|) and taking into account eq. flTUl) one can show that for arbitrary 
(3 we have 



2M 8 
< 9' 



R 



ph 



(12) 



As an illustration of the above presented method for generating exact dark energy star 
solutions we will give a fully explicit dark energy star interior solution generated from 
the well-known interior Schwarzschild solution |28j . The interior Schwarzschild solution 
qualitatively describes the general case of a static, spherically symmetric perfect fluid star 
in general relativity and, therefore we expect that the generated solution should qualitatively 
describe the general case of mixed dark stars. With the generation techniques applied on 
the interior Schwarzschild solution we obtain: 



; m-Vm 2 -d 2 \ 



dr 2 



P 



3VM 2 - D 2 2 M -/f- p2 A 

-g VA/ 2 -D 2 



AttR 3 



1/2 



+ r 2 (d6 2 + sin 2 9d(j) 2 ) 



i _ 2VM 2 -D 2 r 2 
1 /? 3 ' 



^ _ 2VM 2 -D' 2 



1/2 



3 1 



2VM 2 -D 2 
R 



1/2 



2VM 2 -D 2 o 
r 



1/2 



3M 2 m-v^CT a M - VM 2 - D 2 
p = — Vm 2 -d 2 _|_ 3 ^ = ^_ — p ; 



AttR 3 



3D 2 M -/y 2 -f 2 A 

p D = ^ o e Vm 2 -d 2 _|_ 



VM 2 - D 2 
3D 



AttR 3 
D 

VM 2 - D 2 



VM 2 - D 2 



A. 



,(13) 
(14) 

(15) 

(16) 
(17) 
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where 



In 



2VM 2 - D 2 
R 



1/2 



-2 11 



2VM 2 - D 2 



1/2- 



(18) 



Contrary to the pressure, the fluid energy density p and the dark charge density p^, do 
not vanish on the boundary just as in the case of the fluid energy density in the interior 
Schwarzschild solution. We have 



p(R) = 
Pd(R) 



3M 
3D 



M-\/ M 2 -D 2 

2y/M 2 - D 2 \ V^ 1 ^ 1 
R 



AltR* 



1 - 



M-\/ M 2 -D 2 

2VM 2 - D 2 \ 
R 



(19) 



(20) 



Our interior solution is completely regular everywhere for < r < R if the dark star 
radius R satisfies the inequality 



\JM 2 - D 2 4 
R < 9' 



(21) 



On the dark star surface r = R, our interior solution ffl3|) matches continuously the 
exterior solution (El) which in terms of M and D reads 



ds 



ext 



2VM 2 -D 2 \ 



I 1 - 
D 



r°ext 



2VM 2 - D 2 



2VM 2 - D 2 



In 1 



M-\/M 2 -D 2 
\/M 2 -D 2 



2VM 2 - D 2 



dr 2 



2VM 2 -D 2 
r 



+ r 2 {d0 2 + sin 2 



,(22) 



(23) 



This solution reduces to the exterior Schwarzschild solutions in the absence of dark 
energy (i.e. for D — 0). 

We will also consider extremal dark star configurations which can be obtained from (fl3|) 
by taking the limit |D| — > M and keeping R fixed. The described limit gives the following 
interior solution 



dsl 



P 



3- 



^)dt 2 + e~^V~^) [dr 2 + r 2 (d6 2 + sin 2 Odcj) 2 )] , 



47T/? 3 



PD = ±P, 

p = 0, 



(24) 



This extremal interior solution matches continuously the exterior solution which is obtained 
from (1221) in the limit \D\ — > M, namely 

ds 2 e = -e^dt 2 + [dr 2 + r 2 {d9 2 + sin 2 6d<p 2 )] , (25) 

M , s 

= =F — • (26) 
r 

As one can see from ( |24l) the pressure p of the ordinary matter vanishes in the extremal 
limit. How can one explain the existence of such extremal dark stars with zero pressure? 
Generally speaking, the phantom field yields repulsive rather than attractive force. From 
the contracted Bianchi identities 

d r p +(p + p)d r U = p D d r (p, (27) 

describing the equilibrium, one can see that the dark energy term pnd r (p provides additional 
effective pressure which can balance the gravitational forces. 

The extremal configurations in the general case can be investigated on the basis of the 
field equations ([T]). Following |23j one can show that for extremal configurations with zero 
pressure (p = 0) the 3-metric = e~ 2U gij(i,j = 1,2,3) is flat (i.e = 5ij), po = ±p, 
if = ±U and that U satisfies the equation 

A f U = 47rpe- 2l/ , (28) 

where Af is the ordinary flat Laplacian. Solving this equation in vacuum, i.e. for p — 0, 
and taking into account that = 5ijC 2U and tp = ±U we obtain the exterior solution (|25|) . 
Therefore the exterior extremal solution (1251) is the same for all extremal configurations 
with zero ordinary matter pressure . 

Exact extremal interior solutions can be found by specifying the dependance p = p(U). 
For example the extremal interior solution (12^1) is obtained for pe~ 2U = const. Another 
exact and physically well behaved solution is found for pe~ 2U = —C 2 /Att x U where C > 
is a constant. In this case we have to solve the equation 

A f U + C 2 U = 0. (29) 

The solution of this equations which is well behaved at the center r = is 

U = A^, (30) 
r 

with A being an integration constant. In order to match the exterior and interior solution 
we require that the potential U and its radial derivative are continuous on the boundary 
r = R. Imposing these conditions we obtain the following solution 



U 



-M- 



sm 



/ 7r r_\ 
\2 RJ 



(31) 
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Choosing appropriate functions p = p{U) one can construct many other extremal solutions. 

Finally, we give the upper bound of the ratio 2M/R p h for the extremal configurations, 
namely 



2M 



2M 



g -m/r < 2 



(32) 



R 



e 



In conclusion, the exact interior solutions (and the exterior solution) found in the present 
paper could be used in studying, both qualitatively and quantatively, the local astrophysical 
effects related to the existence of dark energy. Some implications of the solutions presented 
in the present paper will be considered in a future publication. There remain some issues 
that need further investigation. One such issue is the stability of the found solutions. 
The fact that the phantom scalar field exhibits more stable behavior than its canonical 
counterpart [T7], [18] and the results in [13] show that we should expect that the non- 
extremal solutions are classically stable when the seed Einstein-perfect-fluid solution is 
stable. This problem will be considered more thoroughly in a future publication. 

The present work could also be extended in the following direcrion. Performing an 
approprate analytical continuatioin of the exterior solution ((6]) we obtain the well-known 
Bronnikov-Ellis phantom scalar (vacuum) wormhole solution [29], [30]. It is interesting 
whether Bronnikov-Ellis branch can also be the external field of a star and what is the 
structure of such a star. We anticipate that some interior solutions countiniously matching 
the exterior Bronnikov-Ellis solution could be obtained by an approprate analytical con- 
tinuatioin of some of the interior solutions constructed in this paper and probbaly some of 
those solutions would have structure similar to that considered in [T5] . 
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